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摘 要 : 赴 线 性 系统 辨识 趾 现 代 辨 识 领域 中 的 一 个 主要 问题 。 在 非 线 性 系统 辨识 中 ， 系 统 常 被 表示 为 一 系 
列 块 连接 。 和 针对 砷 线性 系统 中 的 Hammerstein 模型 ， 本 文 所 出 了 利用 混合 粒 了 和 群 优化 算法 对 非 
线性 系统 模型 进行 辨识 。 该 方法 的 基本 思想 是 将 非 线 性 系统 的 辨识 问题 转化 为 参数 空间 上 的 优化 
问题 ， 然 后 采用 粒子 群 优化 算法 (PSO) 获得 该 优化 问题 的 解 。 为 了 进一步 增强 粒子 群 优 化 算法 
的 辨识 性 能 ， 提 出 利用 一 种 混合 粒子 群 优 化 算法 。 最 后 ， 给 出 仿真 实验 ， 其 结果 验证 了 本 文 给 出 
的 辨识 方法 是 有 效 的 。 
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1 引言 


线性 系统 的 辨识 研究 理论 已 趋 于 成 熟 &3 ， 但 在 现实 生活 中 ， 实 际 系统 几乎 是 非 线性 系 
统 ， 所 以 对 非 线性 系统 的 研究 十 分 必要 加。 然而， 目前 尚 缺 少 描述 各 种 非 线 性 系统 特性 的 
统一 数学 理论 ， 因 而 非 线性 系统 的 辨识 往往 是 和 特定 的 非 线性 系统 模型 描述 相对 应 和 内。 在 实 
际 应 用 中 ， 常 常 许多 非 线 性 系统 可 用 无 记忆 的 非 线性 增益 环节 与 线性 子 系统 的 内 连 来 模型 
1L 45, Hammerstein 模型 是 一 种 非 线性 系统 的 模块 模型 ， 它 是 由 一 个 静态 非 线 性 和 一 个 动态 
线性 模块 串联 构成 ， 大 部 分 实际 非 线 性 系统 可 用 Hammerstein HAKR., 

粒子 群 优化 (Particle Swarm Optimizer, PSO) 算法 是 一 种 进化 计算 技术 ， 最 早 是 由 Kennedy 
fll Eberhart 于 1995 人 年 提出 四， 它 可 以 用 于 解决 非 线性 、 不 可 微 和 多 峰值 的 复杂 问题 。 由 于 其 
思想 简单 、 操 作 易 行 和 解决 问题 的 有 效能 力 而 被 应 用 到 许多 领域 由。 当然 ， 它 也 为 解决 复杂 非 
线性 系统 的 辨识 问题 提出 了 一 条 可 能 的 途径 。 但 在 实际 应 用 当中 ， 也 表现 出 了 一 些 不 尽 人 意 的 
问题 ， 这 些 问 题 中 最 主要 的 是 它 容易 产生 早熟 收 你 、 局 部 寻 优 能 力 较 养 等 &9] 。 

本 文 利 用 在 一 般 PSO 算法 搜索 过 程 中 融入 确定 性 局 部 搜索 算法 和 变异 操作 算法 ， 以 此 构成 
一 种 混合 粒子 群 优化 算法 (Hybrid Particle Swarm Optimization, HPSO) 来 克服 算法 的 早熟 收 
敛 、 局 部 导 优 能 力 较 凑 等 等 缺点 ， 将 该 算法 用 于 对 Hammerstein 模型 的 辨识 ， 仿 真 结果 说 明 其 
辨识 精度 明显 提高 ， 收 到 了 较 满意 的 效果 。 


2 一般 粒子 群 优化 算法 


PSO 是 粒子 在 解 空 间 追 缚 最 优 粒子 进行 搜索 。 它 首先 在 可 行 解 空 间 中 随机 初始 化 一 群 粒 
子 ， 每 个 粒子 都 为 寻 优 问题 的 一 个 可 行 解 ， 并 用 某 个 函数 计算 出 相应 的 适应 值 以 确定 是 个 达到 
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寻 优 目标 。 每 个 粒子 将 在 解 空间 中 运动 ， 并 由 一 个 矢量 决定 其 运动 方向 和 人 位移。 通常 粒 子 将 追 
随 当 前 已 知 的 最 优 位 置 ， 并 经 逐 代 搜 索 ， 最 后 得 到 最 优 解 。 
设 第 ;个 粒子 位 置 为 zi = [xil, Xi2,… Ein] EEK v = [vil, vi ,uia 工 。 它 迄今 搜索 
到 的 个 体 极 值 和 种 群 极 值 分 别 为 p; = [pas Pi2; Pin] T> Pg = [pg1i, Pg2,"… ,Pon] 。 按 追随 当 
前 已 知 的 最 优 位 置 的 原理 ， 粒 子 i 将 按 (1) 式 改变 其 位 移 方向 和 步 长 
Via(t + 1) = w(t)via(t) + eiri (pia(t) — via(t)) + c2r2(pga(t) — zialt)), (1) 
ig(t +1) = zia(t) + via(t + 1). 
式 中 4 = 1,2,… ,n。n 为 解 空间 的 维 数 ; i = 1,2,… m. m 为 种 群 规模 ， 一般 取 为 20 到 40; ri 
和 72 为 [0,1] 之 间 的 随机 数 ，c1 和 ec 为 位 移 变 化 的 限定 因子 ， 通 常 取 为 2; 为 惯性 权重 ， 通 常 
采用 线性 递减 策略 ， 即 
w= (Wini F= Wend) (Tmax 一 t)/Tmax + Wend- (2) 


其 中 Thax 为 最 大 进化 代数 ，wini 为 初始 权 值 ，wena 为 进化 至 最 大 代数 时 的 权 值 ， 为 使 粒子 位 


3 ”混合 粒子 群 优化 算法 


3.1 Hooke-Jeeves 模式 搜索 优化 

Hooke-Jeeves 直接 搜索 法 是 由 两 种 移动 组 成 的 序列 优化 技术 ， 一 种 是 探索 移动 ， 另 外 一 种 
是 模式 移动 ， 从 几何 意义 上 讲 ，Hooke-Jeeves 方 法 是 寻找 具有 较 小 函数 值 的 “山谷 ”， 力 图 使 
迭代 产生 的 序列 沿 “ 山 谷 ” 走 向 并 近 极 小 点 。 算 法 从 初始 点 开始 ， 探 测 移动 依次 沿 %w 个 坐标 轴 
进行 ， 用 以 确定 新 的 基点 和 有 利于 函数 值 下 降 的 方向 。 模 式 移动 沿 两 个 相 邻 基点 连 线 方向 进 
行 ， 试 图 沿 着 “山谷 ”使 函数 值 更 快 的 减 小 。 

3.2 变异 操作 

进化 算法 当中 ， 变 异 算 子 可 增加 种 群 多 样 性 ， 从 而 完成 探测 过 程 。 当 然 ， 变 异 操作 的 引入 
使 PSO 算法 的 收敛 速度 下 降 ， 但 是 解 的 精度 可 明显 提高 。 变 异 过 程 可 描述 为 ， 对 每 个 粒子 赋 
一 个 [0,1 之 间 的 随机 数 ， 若 该 随机 数 小 于 给 定 的 变异 域 值 yw,， 则 粒子 重新 撤 向 搜索 空间 ， 其 
记忆 自身 最 优 位 置 和 速度 向 量 保 持 不 变 。 从 而 增强 了 PSO 算 法 全 局 探测 能 力 。 

3.8 混合 算法 

混合 算法 [19] 实现 的 基本 思想 为 : 设想 在 一 维 的 搜索 空间 ， 有 ma, 2o, ze 和 za 四 个 局 部 最 优 
点 ， 全 局 最 优点 在 z。 如 果 两 个 粒子 a 和 5 在 各 自 的 邻 域内 通过 Hooke-Jeeves 方 法 搜索 到 局 部 
的 最 优 值 ， 然 后 再 利用 PSO 算法 自身 比较 特征 ， 则 全 局 最 优点 很 容易 即 可 找到 。 

在 优化 过 程 中 ， 局 部 搜索 的 嵌入 有 两 种 方式 。 其 一 ,每 隔 给 定 代数 下 后 (比如 了 =1,10 或 100)， 
对 所 有 粒子 实施 局 部 搜索 ;其 二 ， 每 代 都 实施 局 部 搜索 ， 但 不 是 对 所 有 粒子 ， 而 是 设 定 一 个 局 
部 搜索 概率 域 值 jwa, ua € (0,1)。ja 来 控制 实施 局 部 搜索 粒子 的 多 少 。 在 第 二 种 方式 中 ， 还 可 
以 只 对 最 优 粒 子 进行 局 部 搜索 ， 或 对 每 个 粒子 自身 记忆 最 优点 进行 局 部 搜索 ， 若 个 体 自身 记忆 
最 优点 没有 改进 ， 则 无 需 进行 局 部 搜索 。 同 样 ， 变 异 操作 也 有 两 种 方式 ， 车 jxa = 0， 则 粒子 不 
执行 Hooke-Jeeves 操 作 ， 此 时 若 um = 0， 则 HPSO 成 为 PSO。 

HPSO 算 法 描述 如 下 。 

步骤 1 设 定 种 群 m，Hooke-Jeeves 搜 索 概 率 域 值 jJa， 变 异 概 率 域 值 jw,， 最 大 进化 代 
数 Tmsw， 惯 性 权 值 win; 和 wena， 进 化 代数 计数 器 t = 1， 初 始 位 置 (zx1(t), x2(t),…… ,zm(t))， 
初始 速度 (vi (t), vo (t), - - ,Vm(t)); 
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ERA 评价 各 粒子 适应 值 ， 对 第 i 个 粒子 配置 分 布 于 [0,1] 随 机 数 ， 若 其 小 于 ja， 则 按 以 
下 执行 Hooke-Jeeves 搜索 ， 人 奋 则 ， 进 行 步 又 3，; 

步骤 2.1 当前 点 zi(t) 执行 Hooke-Jeeves 搜 索 ， 在 每 个 坐标 方向 (el, e2,-… ,ep) 均 选 定 初 
始 步 长 65， 并 设 定 上 = 1; 

步骤 2.2 计算 所 二 f(z4+D( 引 )。 设 j 21, y = xz 从， 按 步骤 2.3 开 始 探测 移动 ; 

步骤 2.3 ”根据 当前 基点 yO 获得 新 的 基点 ， 以 使 变量 cU? 人 的 扰动 


yO 8e; F= f(y + 6e;) < f — f(y(), 
yit = 4 y se F= flO 5e) < f = f(y?), (3) 
f = f(y0)) < min(f, f). 


这 种 寻找 新 基点 的 过 程 持续 7 = 1,2,-., EISE UP (0) 被 扰动 而 找到 WP+D; 

步骤 2.4 若 WD+D 仍 保持 zt 人 ， 减 少 步 长 5 = B66， 其 中 6 为 缩减 率 ，6 e (0,1)。 设 j= 
1， 到 步骤 2.3。 若 JWD+D 与 zt9( 伯 不 等 ， 则 新 基点 zk+0( = yYD+D， 到 步骤 2.5; 

步骤 2.5 根据 z(t) 和 z+D(t)， 建 立新 模式 移动 方向 s，s = aD (t) -Pe H 
Si y (D, y = z+D(t) + as， 其 中 a 为 步 长 加 速 因子 

步骤 2.6 Wk—kc1,j—1, fu — jwG)， 重 复 步 骤 2.3。 若 jyCD+D) < f(z f 
新 基点 z+D(t) = y(D+D， 到 步骤 2.5。 车 f(y(D+D) > f(a(P (t), Bat = «(9 
小 步 长 6= pi, 令 k= 上 十 1， 到 步骤 2.2; 

步骤 2.7 如 此 继续 下 去 ， 直 到 满足 精度 要 求 ， 即 6 < se。 

步骤 3 比较 粒子 的 适应 值 和 自身 最 优 值 pbest， 如 果 当 前 值 比 pbest 更 优 ， 则 置 pbest 为 当 
前 值 ， 并 设 pbest 的 位 置 为 马 维 空间 中 的 当前 位 置 ; 

步骤 4 比较 粒子 适应 值 与 种 群 最 优 值 gbest， 如 果 当 前 值 比 gbest 更 优 ， 则 置 gbest 为 当前 
粒子 的 矩阵 下 标 和 适应 值 ; 

步骤 5 按 (3) 式 更 新 方向 和 步 长 ， 得 新 种 群 关 (t 十 1)， 检 查 速 度 和 位 置 是 否 超 出 限制 ; 

步骤 6 执行 变异 操作 ， 对 第 i 个 粒子 配置 [0,1] 随机 数 ， 若 其 小 于 jww,， 执 行 变异 操作 ; 

设 t= 二 tt 十 1， 重 复 步 骤 2 到 步骤 6 直到 满足 停止 标准 。 


t))， 得 
nt), I 


4 Hammerstein 模型 辨识 


Hammerstein 模型 由 是 由 一 个 无 记忆 非 线性 增益 环节 和 线性 子 系统 串联 而 成 ， 其 连接 方式 
如 图 1 所 示 。 该 文 考虑 一 类 有 色 噪 声 干 扰 下 的 Hammerstein 模型 的 差分 方程 可 描述 为 


A(q ^ )y(k) = B(q a(k) + C(a- wlk), (4) 
其 中 g-! 为 灌 后 算 子 ，z(k) 为 非 线 性 增益 环节 输出 ，y(%) 为 线性 子 系统 输出 ，w(k) 是 均值 为 
F, JAAP 的 高 斯 白 噪 声 ，C(9 一 ) 为 白化 滤波 器 ， 且 w(k) 与 u(k) 相互 独立 。 且 


n m l 
A(q)-14- aq 5, B) =Y b, C(a))-1-4 Y, cia 
i—l 24=1 i=l 
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无 记忆 非 线性 增 蔓 


图 1: Hammerstein 模型 结构 图 


无 记忆 非 线性 增益 环节 采用 规格 化 以 后 的 pp 阶 指数 多 项 式 有 逼近 ， 即 
p 
x(k) = u(k) + V ^ riu'(k). 
j—2 
且 假设 系统 (4) 的 特征 方程 z? 二 a1z"71 + a2z^7? 十 .… 十 an 一 0 的 特征 根 都 在 单位 圆 以 内 。 
引入 参数 向 量 
0 = (Q1, a2, ica , ân, b1, b2,* -- , bm, C1, €2,* t ,€65T$2,T3, 777 ,7p)T， 
则 辨识 的 目的 就 是 在 给 定 输入 信号 ulk) 和 系统 输出 y(k) 的 情况 下 估计 参数 向 景 9。 
设 参 数 向 量 6 的 估计 值 
6 二 (à1,2,--- ûn, b1, b2, + Ôm, 61,62; zn , Ĉt, fa, f3,- 2s fo). 
则 估计 的 偏差 可 以 用 以 下 准则 函数 外 来 衡量 
J(k) = > 8( o [wk — i) - 9( - 2], (5) 
i=0 


其 中 s 为 辨识 窗口 宽度 ，B(k,i) HRE, 0k) 为 根据 估计 模型 计算 出 的 输出 值 。 
由 于 极 小 化 (5) 式 是 一 个 优化 问题 ， 可 以 采用 以 上 给 出 的 混合 粒子 群 优化 (HPSO) 算法 来 
求 (5) 式 的 极 值 及 其 对 应 的 模型 参数 值 。 


5 ”仿真 示例 


考虑 系统 形式 如 (4) 式 的 Hammerstein 模 型 ， 其 中 喉 声 w(k) 的 均值 为 零 ， 根 方差 o = 
0.01， 输 入 信号 w(k) 为 零 均值 高 斯 白 噪 声 序列 ， 方 差 o? — 1， 取 
A(q 1)-1-—1.5q971 0.797, B(q71) =q! + 0.5q7?, 
C(q |) 21-15g !, x(k) = u(k) + 0.5u?(k) +0.3u3(k) + 0.1u*(k). 
由 以 上 可 知 ， 该 模型 总 共有 8 个 参数 需要 估计 ， 即 9 = (a, a2, bi, b2, c1, T1, T2 T3, T4)] - 
辨识 算法 参数 设置 窗口 s = 50， 对 任意 k 和 i， 权 重 B(k,i) = 1， 采 用 HPSO 对 (5) 式 寻 


优 ， 粒子 数 m = 15， 适 应 值 精度 sy = 0.001， 初 始 和 最 终 惯性 权 值 wis; = 0.9，wena = 0.4, 
信任 参数 cl = ca — 2， 最 大 速度 Vnax = 2，HJ 搜 索 的 概率 域 值 we = 0.05， 变 异 的 概率 域 
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fiu, = 0.05， 初 始 步 长 6 = 0.5， 加 速 因子 a = 1， 缩 减 率 8 = 0.5， 人 允许 精度 e = 0.01, WR 
次 数 为 50 和 所 有 参数 搜索 范围 初始 化 为 [一 2, 2]。 

表 1 列 出 利用 HPSO 算 法 和 PSO 算 法 进行 随机 仿真 试验 10 次 的 统计 结果 。 从 表 1 中 所 列 结 
果 可 以 看 出 ， 利 用 HPSO 算 法 的 辨识 方法 辨识 的 精确 度 明 显 提 高 ， 充 分 显示 了 应 用 HPSO 算法 
的 辨识 方法 的 有 效 性 和 可 行 性 。 且 此 时 应 用 HPSO 算法 辨识 得 到 的 相应 估计 模型 与 真实 模型 的 
阶 跃 响应 如 图 2 所 示 。 


表 1: 10 次 随机 仿真 试验 的 统计 结果 


HPSO 估计 值 PSO 估 计 值 


表示 更 沦 值 
* _ 表示 估计 信 


输出 数值 


miats 


图 2: HPSO 算法 辨识 的 阶 跃 响应 


6 结论 


本 文 利用 了 一 种 混合 粒子 群 优化 算法 对 Hammerstein 模型 进行 辨识 。 混 合 粒子 群 算 法 由 于 
局 部 搜索 能 力 的 增加 ， 克 服 了 早熟 收敛 问题 ， 可 找到 问题 的 最 优 解 。 仿 真 结果 表明 ， 混 合 粒子 
群 优化 算法 对 非 线 性 系统 的 辨识 精度 明显 提高 ， 显 示 了 该 方法 的 可 行 性 和 有 效 性 。 
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Abstract: The nonlinear system identification is a main topic of modern identification. In the non- 
linear system identification, a system is often represented as a series of blocks. According to the 
Hammerstein model, this article proposes a method for the nonlinear system identification by using 
the hybrid particle swarm optimization (HPSO). The basic idea of the method is that the problem of 
nonlinear system identification is changed into an optimization problem in parameter space and the 
particle swarm optimization (PSO) is then adopted to solve the optimization problem. In order to 
enhance the performance of the PSO identification, an HPSO is also presented. Finally, simulation 
results show that proposed methods are effective. 
Keywords: system identification; hybrid; particle swarm optimization; Hammerstein model 
Received: 09 May 2007. Accepted: 05 Mar 2008. 
Foundation item: The National Natural Science Foundation of China (60874033); the Fund for the 


Doctoral Program of the Education Department of China (20060700007); the Project Item of Technological 
Innovation of Xi'an University of Technology (108-210813); the Academic Fund of Xi'an Jiaotong University 


(0109-08140015). 


